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Abstract 

The free energy and local height probabilities of the dilute A models with broken 
7L2 symmetry are calculated analytically using inversion and corner transfer matrix 
methods. These models possess four critical branches. The first two branches provide 
new realisations of the unitary minimal series and the other two branches give a direct 
product of this series with an Ising model. We identify the integrable perturbations 
which move the dilute A models away from the critical limit. Generalised order pa- 
rameters are defined and their critical exponents extracted. The associated conformal 
weights are found to occur on the diagonal of the relevant Kac table. In an appropri- 
ate regime the dilute A3 model lies in the universality class of the Ising model in a 
magnetic field. In this case we obtain the magnetic exponent 5 = 15 directly, without 
the use of scaling relations. 
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1 Introduction 



In the last decade many infinite hierarchies of exactly solvable models have been found. Of 
foremost importance among these models are the restricted solid-on-solid (RSOS) models of 
Andrews, Baxter and Forrester (ABF) [|[]. In these models each site of the lattice carries a 
height variable, restricted to the values 1, . . . , h — 1 with h = 4,5, . . ., subject to the rule 
that heights on neighbouring lattice sites differ by ±1. If the allowed heights are represented 
by the nodes in the following diagram. 



h-l 



the adjacency rule requires that neighbouring lattice sites take values that are adjacent on 
the diagram. 

Andrews, Baxter and Forrester considered four different regimes, labelled I-IV. It was 
pointed out by Huse that the critical line separating regimes III and IV realises the 
complete unitary minimal series of conformal field theory. This series is described by a 
central charge 

and a set of conformal weights, given by the Kac formula 

^- = ""'^M/^-'lf'' l<r<h-2, l<s<h-l. (1.2) 

The corresponding modular invariant partition function is 

^ = ^EEIxS(^)l' (1-3) 

^ r=l s=l 

where q is the modular parameter and the xi'^] ^ire the Virasoro characters given by 

a''^'— c/24 oo 

^ ^ g ^^h{h-l)f + [hr-{h-l)s]j_^h{h-l)f + [hr+{h-l)s]j+rs)^ (^^^^^^ 

Q{q) j=-oo 

withg(g) = n^=i(i-g"). 

By giving a loop or polygon interpretation to the critical ABF models, Pasquier ^ 
extended these models to arbitrary adjacency graphs. Demanding that these new models be 
critical restricts the graphs to the Dynkin diagrams of the classical and affine simply-laced 
Lie algebras shown in Fig. 0. 

Recently a new construction of solvable RSOS models was found 0, |^, |^. Basically, the 
method is an extension of the work of Pasquier, and related work of Owczarek and Baxter 
1^, to more general loop models. Application to the 0(n) model [|10], which is closely related 
to the Izergin-Korepin model has led to a new family of critical RSOS models labelled 
by Dynkin diagrams. The same models were found independently by Roche [|I2|. 



In the approach of Pasquier, the polygons, which are interpreted as domain walls separat- 
ing regions of different height, densely cover the edges of the dual lattice. As a consequence, 
heights on adjacent sites are always different. In the new RSOS models, two neighbouring 
sites of the lattice either have the same or different height, so that the domain walls occupy 



some but not all edges of the dual lattice. Therefore it is natural, following [|T^, to term 
these new models dilute A-D-E models. 

Each member of the dilute A^ hierarchy possesses four distinct critical branches. The 
central charge is given by 

■j^jj^ — branch 1 and 2 



3 6 



2 h{h-l] 



(1.5) 
branch 3 and 4 



where 



^ _ f L + 2 branch 1 and 3 ,^ 

~ \ L + l branch 2 and 4. ^ ' ^ 

The first two branches give new realisations of the unitary minimal series with the modular 
invariant partition functions (|1.3| ). The other two branches appear to be a direct product of 
this same series and an Ising model, with modular invariant partition functions 

^ = I E E EE Ixit(g) x^M'- (1-7) 

^ r'=l s'=l r=l s=l 

As reported in 0], the models related to the Al Dynkin diagrams admit an off-critical 
extension. A remarkable feature of these off-critical models is that, for odd values of L, 
they break the Z2 symmetry of the underlying Dynkin diagram. The simplest of these 
symmetry breaking models belongs to the universality class of the Ising model. This allows 
the calculation of the magnetic exponent 6 = 15 without the use of scaling relations. 

This paper is devoted to the investigation of the models of the dilute Al hierarchy. First 
we briefly describe the whole family of dilute A-D-E models. Then, in section ^, we define 
the off-critical A^ model and in section ^ we calculate its free energy. From this we extract 
the critical exponent a when L is even and 6 when L is odd. The main body of the paper 
is concerned with the calculation of the order parameters of the dilute A models for odd 
values of L. In section |^ we compute the local height probabilities and in the subsequent 
section we use these results to evaluate generalised order parameters. We also extract the 
set of associated critical exponents 6k and derive the corresponding conformal weights. In 
section |^ we discuss the phase diagram, concentrating on L = 3, and in section ^ we collect 
results concerning the Ising model in a field. Finally, we summarise and discuss our main 
results. 

The results for the order parameters when L is even will be presented in a future pub- 
lication. Likewise, results for the critical models related to the other adjacency diagrams, 
among which is a solvable tricritical Potts model |T^, will be reported elsewhere. 



2 The dilute A-D-E models 



In this section we define the family of dilute A-D-E models. Although we restrict the 
description to the square lattice, they can be defined on any planar lattice. 

Consider an arbitrary connected graph Q consisting of L nodes and a number of bonds 
connecting distinct nodes. Label the nodes by an integer height a G {1, . . . ,L}. Nodes a 
and b are called adjacent on Q if they are connected via a single bond. Such a graph is 
conveniently represented by an adjacency matrix A with elements 

. _ J 1 if a and b adjacent , . 

"jo otherwise. ^ ' 

Let n denote the largest eigenvalue of A and S the Perron- Frobenius vector, i.e., AS = nS. 

With these ingredients we define an RSOS model on the square lattice C as follows. 
Each site of C can take one of L different heights. The Boltzmann weight of a configuration 
is non-zero only if all pairs of neighbouring sites carry heights which are either equal or 
adjacent on Q. The weight of an elementary face of the RSOS model is given by 

Pl^a,b,c,d + P2^a,b,c^a,d + P3^a,c,(i^a,b + ( TT ) Pih,c,d^a,b 




s, 



+ [^J P5Sa,b,dAa,c + Pfs5a,b5c,dAa,c + P75a45h,cAa,b (2.2) 



+ P?,5a,cAa,b^a4 + (^^gT^ j P9h,d^a,bAh,c 

where Sa is the a-th entry of S and a, 6, c and d can take any of the L heights of the graph 
Q. The generalised Kronecker 5 is defined as = 11^12 ^h,ij- parametrise n by 

n = -2cos4A (2.3) 

then pi, . . . , pq are given byQ 

: (sin 2A sin 3 A + sin u sin(3A — u))/ (sin 2A sin 3 A) 

-- sin(3A — u)/ sin3A 

: sinw/ sin 3 A 

: sinw sin(3A — M)/(sin 2A sin3A) (2.4) 

: sin(2A — u) sin(3A — u)/ (sin 2A sin 3A) 

: — sin-u sin(A — 'u)/(sin2A sin3A). 

We note that the weights p4 = ps and the weights pe = Pi are determined only up to a sign. 
For any graph ^, the RSOS model defined by (|2.2| )- (|2.4| ) satisfies the Yang-Baxter equation 





Pi 


P2 


= P3 


P4 


= P5 


P6 


= P7 




P8 




P9 



H. In fact, in i, |, i it was shown that all models defined above have the same partition 



function as the 0{n) model [p!0|j . 

^We note that we have changed the variable A of reference H to ^tt 



From equations ( p.3|) and ( p.4|) it follows that there are four different branches that yield 
the same values of n. Using the periodicity of the weights, we can restrict A to the interval 



D — i 



branch 1 
branch 2 
branch 3 
branch 4 



< M < 3A 
< M < 3A 
-vr + 3A < M < 
-vr + 3A < M < 



ivr < A < ivr 
ivr < A < ivr 



W < A < iyr 



(2.5) 



ivr < A < ivr. 



These four branches correspond to (part of) the four branches defined in for the 0(n 



model. For n > 2 the value of A — ;|7r must be chosen imaginary, and the weights become 
complex, unlike the ordinary A-D-E models. As was pointed out in 0, 0, the only adjacency 
graphs that have largest eigenvalue n <2 are the Dynkin diagrams of the classical {n < 2) 
and affine [n = 2) simply-laced Lie algebras shown in Fig. |^. For the classical case the 
respective values of n are listed in Table |. The corresponding Perron-Frobenius vector can 
be found in [ffl. 



Classical 



1 2 3 



1 2 3 



2 3 4 5 



1 



2 3 4 5 6 



L-1 



2 3 4 5 6 7 




1 2 3 




^0) 



3 4 



2 3 4 5 



L-1 



L-1 



2 3 4 5 6 7 
9 T 



2 3 4 5 6 7 8 



Figure 1: Dynkin diagrams of the simply- laced Lie algebras. 

From the equivalence with the 0(n) model, the central charge of the dilute A-D-E models 
is known |TB|. The values on the four branches are listed in Table | for the classical algebras. 
For the affine algebras we have X = jn and c = 1 or |. 



3 The off-critical A model 



The dilute A^ model ( p.2|) admits an extension away from criticality while remaining solvable. 
In terms of the theta functions of appendix |^, suppressing the dependence on the nome p, 



algebra 
A. 



E7 



n 



2cos^ 



2 cos 



2L-2 



2cos^ 



2cosf| 



2cos^ 



4 \ + L+l 



^ fl T ^ 



2L-2 



- (1 T - 

4 V ^ 12 



4 \ 18 



4 V ^ 30 



1 



1 



1 



1 



(L+1)(L+1±1) 
6 

(L+l){L+lTl) 



(2L-2)(2L-2±1) 

3 6 

2 (2L-2)(2L-2q:l) 



12(12±1) 
3 _ 6 
2 12(12=Fl) 



18(18±1) 

3 6 

2 18(18^1) 



30(30=Fl) 
3 _ 6 
2 30(30ilil) 



branch 

1,2 
4,3 



1,2 
4,3 

1,2 
4,3 

1,2 
4,3 

1,2 
4,3 



Table I: Central charge of the dilute A-D-E models. 
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Boltzmann weights of the off-critical A^, model are given by 

- u)'di{3\ + u) 



a a 
a a 



W\ 



a ± 1 a 
a a 

a a 
a ± 1 a 

a a ± 1 
a a ± 1 



^9i(6A)^?i(3A) 

'S{a + 1) i^i{2aX - 5A) ^ S{a - 1) ^94(2aA + 5A)\ Mu)A{5X - u) 



S{a) M'^aX + X) S{a) M^aX ~ X) J ^9i(6A)i?i(3A) 
i3i{3X-u)di{±2aX + X-u) 



W 



W\ 



a \ _ 

a a±l J 

a a ± 1 \ 
a a I 

'a±l a±l 
a a 



^9i(3A)i94(±2aA + A) 
S{a ± 1) \ ^i{u)^4{±2aX - 2X + u) 



?9i(3A)t94(±2aA + A) 



'^94(±2aA + 3A)^94(±2aA - A)^^^ ^9i(m)?9i(3A - u 



^l{±2aX + A) 



^9i(2A)?9i(3A) 



W\ 



w 



a±l 



a±l 



a =F 1 



a =F 1 



^i(2A-m)^i(3A-m) 
^9i(2A)?9i(3A) 



5(a - l)S{a + l)y/^ ^9i(m)7?i(A - u) 



S^(a) 



i9i(2A)^9i(3A) 



a 

a±l 



a ± 1 \ _ ^9i(3A - u)di{±AaX + 2A + m) S{a ± 1) ■di{u)di{±Aa\ -\ + u) 
a )~ ^9i(3A)^9i(±4aA + 2A) ^ S{a) ^9i(3A)^9i(±4aA + 2A) 

_ ?9i(3A + u)M±4:a\ - 4A + m) 
" ^9i(3A)^9i(±4aA -4A) 

'5(aT l)7?i(4A) 'i94(±2aA - 5A)\ ^9i(M)^9i(±4aA - A + m) 



+ 



S{a) ^i{2X) i?4(±2aA + A) / ^9i(3A)?9i(±4aA - 4A) 



S{a) 



^i(4aA) 
^94(2aA) ■ 



We note that in the critical hmit, p ^ 0, the crossing factors S{a) reduce to the entries of 
the Perron- Frobenius vector Sa'- 



hm 



S{a) Sa 



P-o S{b) Sb ■ 

We also note that, for later convenience, we have relabelled the states of the model a 
L + 1 — a in comparison with those of our earlier definition of the model in H . 



(3.2) 



The Boltzmann weights ( |3.1| ) satisfy the following initial condition and crossing symme- 



try: 



W 



w 



d c 
a h 



d c 
a b 



3X — u 



( S{a)S{c) y/' 
\S{b)S{d)) 



W 



c b 
d a 



u 



and an inversion relation of the form 



u W 



d c 
g b 



-u 



(3.3) 
(3.4) 

(3.5) 



In equation (|2.5| ) four different critical branches were defined. This yields eight regimes 
for the off-critical model 



regime l"*" 
regime 1~ 



reg: 
reg: 

reg: 
reg: 



ime 2"*" 
ime 2~ 

ime 3"'" 
ime 3~ 



regime 4"*" 
regime 4~ 



< p < 1 
-1 < p < 

< p < 1 
-1 < p < 

< p < 1 
-1 < p < 

< p < 1 
-1 < p < 



< M < 3A 



< u < 3A 



3A - vr < M < 



A 



(l + ITT 



(3.6) 



1 

L+l 



3A 



TT 



<n<0 A = f(l-^). 



For regimes 2^ and 3^ we exclude the L = 2 case because the model becomes singular. 



At criticality all A^, models satisfy the ^2 symmetry of the Dynkin diagram, but the 
off-critical models, for odd values of L, break this symmetry: 




^ W 




L odd. 



(3.7) 



For L = 3, if we make the identification {1,2,3} = {+,0, — }, the model can be viewed as 
a spin-1 Ising model. For p the up-down symmetry of this model is broken. We can 
therefore regard the nome p of the ■(^-functions as a magnetic field. This is in contrast with 
the usual role of p as a temperature-like variable (see, e.g., ^^). Also for larger, odd values 
of L we refer to p as the magnetic field, even though it is not, in general, the leading magnetic 
operator. For odd L the -|- and — regimes are equivalent because negating the magnetic field 
p merely has the effect of relabelling the heights a L + 1 — a. 

The dilute Al model defined above is closely related to the Izergin-Korepin (or A2^'') SOS 

with 2aX replaced 



model [|T^]. In the latter model the weights are given by equation (|3]^ 
by 2aX + wq, where Wq is an arbitrary constant, and a G 2. If we set 



A 



k7T 



k e {1,...L,L + 2,...,2L + 1} 



(3. 



4(L + 1) 
Wo = 

we obtain A2 RSOS models based on the A^ algebra. Only when k = L or L + 2 are these 
models physical. Other choices of k correspond to eigenvalues of the adjacency matrix of 
which are not the largest. For all odd values of k the models break the ^2 symmetry of the 
underlying Dynkin diagram. 

Another way to restrict the Ag SOS model is given by 



A 

Wo 
a 



kn 



2{L + 1) 
= i i In p 

e {i,...,L}. 



k E {1,...L} 



(3.9) 



This possibility has been studied by Kuniba in a more general way in his study of A^^-* 
RSOS models fl^. However, this way of restricting the A^2^ SOS model does not give rise 



to symmetry breaking models. 



4 The free energy 

We calculate the free energy or, equivalently, the partition function per site k of the dilute A 
model by the inversion relation method . Because the inversion relation ( ^.5] ) is quadratic 
in p we can restrict ourselves to the regimes with positive p, and parametrise p = exp(— e). 



4.1 Regimes 1+ and 2+ 

In regimes l"*" and 2"*" we assume that n{u) is analytic in the strip < Re(n) < 3A, and may 
be analytically extended just beyond these boundaries. The inversion relation ( ^.51 ) implies 

= wmm ^ ^ 



while the crossing symmetry ( p.4|) translates to 

k{u) = k{3X-u). (4.2) 

We make the conjugate modulus transformation ( p.8| ) and a Laurent expansion of In k,{u) 
in powers of exp(— 27rM/e). Matching coefficients in equations ([4.1|) and (|4.2|) we obtain for 
the free energy 

1 ( ) 2 cosh(5A — 7r)7rfc/e coshvrAfc/e sinhTT-ufc/e sinh(3A — M)7rfc/e [43) 

k sinhvr^/c/e cosh37rAA;/e 

We now take the p —>■ limit in the above expression to obtain the leading critical singularity. 
Using the Poisson summation formula, we find 

lnfi:sing~p"/'^ (4.4) 
If we compare this with 

luKsing-p^"" or In Ksing ~ P^^^'''^ (4.5) 

we find for the critical exponents a and 6 

a = L even f a = |g| L even 

regime 1"^ : \ regime 2"*" : < (4.6) 

5 = ^ L odd [s= ^ L odd. 



When L = 2 in regime 1+, expression ( [4.4|) for the critical singularity has to be multiplied 
by In p. 

4.2 Regimes 3+ and 4+ 

In regimes 3"^ and 4"*" the appropriate analyticity strip is — tt + 3A < Re(u) < 0, and the 
crossing symmetry becomes 

k{u) = k{3X-7t -u). (4.7) 

Performing the same steps as before we obtain 

, , „ cosh(5A — 7r)7rA;/e cosh7rAA;/e sinhTTufc/e sinh(7r — 3A + n)7rA;/e , , ^, 
\nK(u) = -2 > ^ — -r— — —, , , — —. (4.8) 

For the dominant singularity we find, apart from some exceptions we list below. 



Comparing this with ( [4.5|) , we find for the critical exponents a and 5 

2(L-2) T 

■— L even 



regime 3"^ : < 



a 



2(L+4) T 

— L even 

L odd 



L-2 
L-2 



a 



regime 4'*' : 



L+4 



(4.10) 



3(L+2) 



L+4 
3L 



L odd. 



For the cases L = 5 and 8 in regime 3"^ and L = 2 in regime 4+ equation ( [4.9[ ) has to 
be multiphed by Inp. When L = 3,4, 6 and 14 in regime 3"'" and L = 8 in regime 4"*" the 
partition function per site is regular. 



5 Local height probabilities 

In this section, which forms the main part of our paper, we calculate the local height proba- 
bilities of the dilute A model for odd values of L. Since negating the nome p is nothing but 
a reversal of the magnetic field, we can restrict ourselves to the four '+' regimes. (Recall 
that for odd L the Boltzmann weights are symmetric under the transformation p — > —p, 
a ^ L + 1 — a.) 

5.1 Groundstate configurations 

First we describe the groundstate configurations in each of the four regimes. Below we 
depict the set of groundstates by decoration of the adjacency diagram. Comparing the 
various weights in (|3.1|) in the ordered limit, see appendix 0, we find that the following 



types of groundstates occur: 

1. Completely fiat or ferromagnetic configurations. If the height of this configuration 
is a we will denote this state by a solid circle on the adjacency graph Q at node a. 
Conversely, fiat configurations that do not yield a groundstate will be denoted by an 
open circle on Q. 

2. Antiferromagnetic configurations. One sublattice has height a and the other sublattice 
height 6 = a ± 1. This state, together with the state where the heights on the two 
sublattices are interchanged, is denoted by a double bond on Q between the nodes a 
and b. 

If we define the variable I to be 

2 [-^J + 1 regimes 1+ and 4+ 



(5.1) 

^ 2 [^J regimes 2+ and 3+ 



where [ J denotes the integer part, the groundstates of the four regimes are as shown in 
Fig. 0. Thus the total number of groundstates is |(L + 1), |(L — 1), |(L + 1) and |(L — 1) 
for the regimes 1"*", . . . , 4+, respectively. 



regime 1+ • — o — • — -O — • — • — C- • — O — • — O 

12 3 I l+l L-1 L 



regime 2 • — o — • — O- • — o — O — •— • — o — • — O 

12 3 I l+\ L-1 L 



regimes"^ C i • O — • =o — • • — — » — en 



12 3 Z Z+1 L-1 L 



12 3 I l+\ L-1 L 

Figure 2: Groundstate configurations for tlie four different regimes. 

5.2 Local height probabihties 

Tfie local lieigtit probability P^^{a) is tlie probability that a given site of the lattice has 
height a given that the model is in the phase indexed by h and c. In the ferromagnetic 
phases c = h and in the antiferromagnetic phases c = 6 ± 1. 

The technique of corner transfer matrices (CTMs) which is used to calculate P'"^{a) is 
well-known and the details of the method are given elsewhere (see, e.g., 0, IH). Because 



the weights (3J.) satisfy the Yang-Baxter equation and possess the crossing symmetry 



the local height probability in regimes l"*" and 2"^ can be written as 

P'^(a) = hm ^'"^'^^^(^)-^^^^(^) (5 2) 

Here the variable q is related to p = exp(— e) by 

q = e-l2'^^/^ (5.3) 
The one-dimensional configuration sums X^'^{q) are given by 

J^o-i cr,„+i cr,„+2 _ ^ gZ]jllJ-'^('^J'°'j + l'°'^+2)_ {5-4) 

The weight function H is determined from the Boltzmann weights in the ordered limit 
{p — > l,u/e fixed). Their behaviour is 



' a b I gf,g^ 



where 

ga = e-2^"'^'/^ (5.6) 

The values of the function H{d, a, b) for regimes l"*" and 2"*" are listed in appendix 0. 

From ( |5.4| ) it follows directly that the one-dimensional configuration sums satisfy the 
following recurrence relation: 

X^'^iq) = g"^^(^-i'^''=)X^i-ii''(g) + g'"^(^'^'^)X^i\(g) + g™^(*+i'*'^)X^iY'(?)- (5-7) 



The task is to solve this relation, with initial condition 



or, equivalently, 

and with the conditions 



aL+lL 



(g)=0 



(5i 



(5.9) 



(5.10) 



which confine the heights in the recursion relation (|5.7|) to the set {1, . . . , L}. 

In regimes 3+ and 4"*", where the spectral parameter u is negative, the local height prob- 
ability is given by 



lim 



g°'^/"^(a)X°'"^(g) 
Ef=i g"'^/^^(a)X«'-=(g) 



where q is now 

q ^ g-47r(7r-3A)/£_ 



(5.11) 



(5.12) 



As for the regimes 1^ and 2"*", the one-dimensional configuration sums are expressed in H 
by ( |5.4| ). To determine if, we now need to know the powers of e"*"^™/^ in the ordered limit 
of the weights. Relative to regimes 1^ and 2"*", after the appropriate gauge factors g are 



separated, and contributions which cancel in the ratio in equation ( [5 .111 ) are discarded, H is 
effectively replaced —H. We thus need only take g — > 1/g in the regime l"*" (2"'') solution of 
( p.7| ) to obtain the form of the solution in regime 4'^ (3^). Of course, q must still be given 
its proper meaning from ( |5.12|) . 



5.3 Solution of the recurrence relation 

In this section we derive the solution of the recursion relation (|5.7| ) for each of the four 
regimes. In fact, as remarked above, the solution for regimes 3"'" and 4"'' is deduced from the 
solution for regimes l"*" and 2"*". 



Gaussian multinomials 



3efore we present the solution we need some preliminaries on the Gaussian multinomials 



defined by 



m 




m 


_k,l_ 




k,l_ 



{(i)k{q)i{q)m-k-i 



(5.13) 



where 



(5.14) 



k=l 



In the limit g — 1 the Gaussian multinomials reduce to ordinary multinomials 



lim 



m 
kj 



m 



k\l\{m-k-l)\ \k,lj' ^^'^^^ 

The following identities, which prove useful in solving the recurrence relation, may be derived 
by straightforward manipulation of the definition 



m 
k,l 



m — 1 
k,l 

m — 1 
k-l, I 

m — 1 
k-l, I 



+ q 



m—k 



+ q 



m—l 



m — 1 
k-l,l_ 




k-l 


m 
k,l 


m—l 
_k,l - 1 


+ q' 


m — l 
_ k,l _ 



k+l 



m — l 




m—l 


_ k,l _ 


+ q' 


k,l-l 



(5.16) 
(5.17) 
(5.18) 



Regime 1 



As a first step towards the solution of equation ( |5.7| ) we consider the special limit q = 1- 
In this limit, setting X^^^l) = X^", the recurrence relation reduces to the following com- 
binatorial problem 

(5.19) 



X. 



abc 



x^ab— 16 1 -yahb _i_ -y 



ab+lb 
m—l 



The fundamental solution of this equation is 



X. 



abc 




(5.20) 



fc=0 



In order to satisfy the initial and confining conditions (|5.9|) and (|5.10| ) we must take linear 
combinations of this solution in the following way: 



X. 



abc 



E 



m 

k,k + 2{L + l)j- 



m 



k,k + 2{L + l)j + a+bj 



(5.21) 



Another piece of information can be gained by considering the m — >■ oo limit. It was 
noted in [20 1 that, in this limit, the one-dimensional configuration sums for the ABF models 
in regimes III and IV become precisely the characters of the related Virasoro algebra ( |1.4| ). 

Expecting similar m ^ oo behaviour for our configuration sums, we generated large 
polynomials on the computer and multiplied these by Q{q). The resulting polynomials, 
being very sparse, were then easily identified as Virasoro characters, with r = a and s = b 
for b < I and s = b + 1 for b > I, ot as linear combinations of two characters. 

Guided by these two limiting cases, replacing h in ( |1.4|) by L + 2 according to (|1.6| ), 
we make the following Ansatz for the configuration sums, with b < I, which yield a single 
character: 



X. 



abc 



q 



7 m 



E 



q 



{L+2){L+l)p + [{L+2)a-{L+l)b]j+a{j,k) 



m 



k, k + 2{L+l)i + a-b 



_^(L+2)(L+l)p + [(L+2)a+{L+l)b\j+ab+P{j,k) 



m 



k,k + 2{L + l)j + a + b 



(5.22) 



For b > I we make a similar Ansatz. The unknown functions a and /3 are assumed to be 
quadratic in their arguments j and k. They depend imphcitly on a, b and c, as does the 
parameter 7. For the configuration sums that, in the m — >■ 00 hmit, yield linear combinations 
of characters, we make the appropriate linear combinations of the above Ansatz. 

Fitting the Ansatz for small values of m with the correct polynomial, the coefficients in 
the functions a{j,k) P{j,k) and 7 can be determined. Before we present the solution, we 
define the following auxiliary function: 



Kib) = q 



{a—s)a/2 



E 



j,k= — OQ 

{L+2){L+l)p + [{L+2)a-{L+l)s\j+k[k+2{L+l)j+a-b 



m 

fc,fc + 2(L + l)j- 



(5.23) 



_^(L+2){L+l)p + [{L+2)a+(L+l)s\j+as+k[k+2(L+l)j+a+b] 



m 



k,k + 2{L+l)i + a + b\ J 

where we have suppressed the a dependence. This function has the following elementary 
properties 



(5.24) 



It is also convenient to define the four sets 



si = {1,3, ...,/} 

52 = {/ + 1,/ + 3,...,L-1} 

53 = {2,4,...,/-!} (5.25) 

54 = {/ + 2,/ + 4,...,L} 

where I is given by ( [5.1| ). With these definitions the solution of the recurrence relation reads 

FM 
q¥F:^\b) 
q^'F^^t\b) 

Ft{b) + {l-q^)q'F'J:l{b+i) 



X: 



abb-l 
m 



q 



mH{b,b,b-l) ^ 



bes,\{l} 
6 e S2 u {L+l} 
6 e S3 
b & S4 



'Y'a.bb „mH{b,b,b) ^ 



X, 



a 6 6+1 



q 



,m_ff(6,6,6+l) 



X < 



FM 
q¥Ft^\b) 

q^'F'^\b) + {l-q^)q-'^^Ftl{ 
F^(6) + (l-g-)g^F^t?(6+l) 

Ftih) 

q^'F:^\b) 

q^'F:^+'{b) + {l-q"^)q~^'Ftl( 
Fiib) 



b e si 
be S2 
be S3 
b e S4 

be si 
be S2 

bessU {0} 
bes,\{L}. 



(5.26) 



We note that in this solution we have included the terms ^ and X^^"*"^ ^ which, according 
to the confining condition (^.10|) , should be identically zero. Using the simple relations ( p.24| ) 
it follows directly that this is indeed the case. The advantage of including these terms in 
the above way is that it enables us to prove the solution without treating the boundary 
separately. From the definition ( [5.23| ) it also follows immediately that the initial condition 
( ^.9|) is satisfied. 

Proving that ( |5.26| ) is indeed the solution of the recurrence relation is now straightforward 
but rather tedious. Inserting ( |5.26| ) into the recursion relation, using the explicit form of the 
function H as listed in equation ( [A. 4] ), we find that only the following four relations need 
hold: 



^m—l 



q-^+^F'^\{h-l) + q 



6+1/ 



l-q 



m—l 



FM-Fi.,ib-1) 



m— 1^ 



76+2 
m-2 



(6+1) 



be S3 



b G S4. 



(5.27) 



The restrictions on b make these special cases of more general expressions, which we will 
show to be true for all values of b. If we widen their applicability in this way, the relations can 
in fact be combined to give simpler ones. In turn, these relations can be further simplified by 
requiring that they hold term-by-term in j. The resulting pair of relations are much stronger 
requirements than the original set, but still they hold. Defining the function 



'^"/s+k(k+"f—b) 



k=0 



m 

k,k + 'y 



i7S+fc(fc+7+6) 



m 

k, k + 'J + b 



6+1 fb-2 



where we have set 2(L + l)j + a = 7, these two relations are 
flib) - flib-l) = 

flip) - fl_,{b) = q-^ 



-6+1/ 



^■m—l\ fb—2 
)Jm-2 



(6-2) 



(5.28) 

(5.29) 
(5.30) 



The function consists of two polynomials. One, the first term within the curly braces, 
has only positive coefficients and the other, the second term, has only negative coefficients. If 
we demand that the equations ( |5.29| ) and ( p.3(J| ) be satisfied for the 'positive' and 'negative' 
polynomials independently, and we set 7 =F 6 = n, respectively, this yields 



k=0 



k{k+n) 



m 
k,k + n 



k=0 



fc(fc+n+l) 



m 

k, k + n + 1 



fc=0 



fc(fc+?l-l) 



m — l 

k, k + n — 1 



,m+n+l 



E^ 

fc=0 



fc(fc+n+2) 



m — l 

k,k + n + 2 



(5.31) 



for equation (|5.29| ), and 



A:=0 



k{k+n) 



m 
k,k + n 



k=0 



k{k+n) 



m — 1 

k,k + n 



^m—l 



E^ 

k=0 



k{k+n+l) 



m — 1 

k,k + n + 1 



qk{k+n-l) 

k=0 



m — 1 
k, k + n — 1 



(5.32) 



+g'"+"(l - g^^-i) ^ g' 



A:(/c+7i+2) 



fc=0 



m — 2 

/c, /c + n + 2 



for equation ( p.30[ ). The proof of these final two equations is now elementary. Equation 
( ^.31|) follows immediately from ( p.l6| ) and (|5.17| ) if we set / = + ra and / = k + n + 1, 
respectively. Equation ( p.32| ) follows from ( p.l6| ) with I = k + n. 

Regime 2+ 

Finding and proving the solution of the recurrence relation for regime 2+ proceeds along 
similar lines as in regime 1"^. To remove any confusion we will denote the one- dimensional 
configuration sums in regime 2+ by Y^^'^. First of all, the q = 1 limit ( p.21j ) is still valid. 
Generating large polynomials reveals that in this case the configuration sums yield Virasoro 
characters with s = a and r = b for b < I ~ 1, and r = b — 1 for 6 > / + 2, or linear 
combinations of two characters. 

We replace h in ( |1.4| ) by L + 1 according to ( |1.6D and define the auxiliary function 



q 



(L+l)Lf-[{L+l)r~La]j+k[k+2{L+l)j+a-b] 



m 



k,k + 2{L + l)j + a-b 



(5.33) 



{L+l)Lj^ + [(L+l)r+La]j+ar+k[k+2(L+l)j+a+b] 



-q 

which satisfies the following simple relations 

Gl^ib) = -G-„^{-b) 

Gt^^\L + l+b) = -g-(^+i)^G^-'-(L + 1-6) 

GliO) = (l-g™)G^_i(l). 

Again we define four sets 

h = {1,3,...,/-!} 

t2 = {/ + 2,/ + 3,...,L- 1} 

ts = {2,4,. ..,0 

U = {/ + 1,/ + 3,...,L}. 



m 



k,k + 2{L + l)j + a + b 



(5.34) 



(5.35) 



The solution of the recurrence relation then reads 



abb-l 



mH(b,b,b~l) ^ 
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{Glib) 
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6 G t2 U {L + l} 
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b et4 
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bet2 
bets 
b et4 

beti 
bet2 
bets 

bet,\ {L}. 



(5.36) 



As in regime l"*" we have extended the solution to include the term Y^^~^^^, which, using 
( p.34| ), indeed yields zero. However, in contrast to regime l"*", we can no longer extend the 
above solution to the term Y^'^^ which according to (|5.36|) would not be zero. Therefore, 
in proving the recurrence relation we have to treat cases that involve this boundary term 
separately. This 'irregularity' is however compensated by the exceptional terms Y^"^^ of 
( p.36| ) and iJ(l,2, 1) of ( |A.5| ). As a result, inserting the above solution into the recurrence 
relation, using ( p.34| ), we again obtain only four relations that should hold: 

Gl{b)-Gl^,{b) 



^m—l 



Gt_,{b+1) + q-'-^'G'^\{b-l) + g^+i(l-g--i)G^„t22(&+2)l b e t. 
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T,m— 1 



G'^\{b-l) + g^+iG^t.\(6+l) + q-'+\l-q^^-')G'^l,{b-2) 



Gt\b) - G''^l,{b-l) 



-,m— 1 



G'-^,{b) + g^+iG^t\(6+l) + g^(l-g--i)G:;t^(6+l) 



Gl{b) - Gl_,{b+l) 



m—l 



Gl^,{b)+q~''+'G'^',ib-l)+q-\l-q-^-')G':^'2ib-l) 



be to 



bets 



betA 



(5.37) 



Again these equations hold for all values of b and we drop the restrictions on b. Doing 
so, the four equations can be combined to give simpler equations. As before, these are true 
term- by-term in j. Setting 2(L + l)j + a = 7, we find 



fn.ib) - f^{b+l) = q^' 



fm{b) - fl^,{b) = q 



m—l 



q^'frn-iib-l)-q''^'CMb+2) 
fm-iib+l) + q-'^' ft;\{b-l) + g^+^(l-g-^)/r_y6+2) 



(5.38) 



with defined in ( ^.28 ). Finally, making the same splitting into 'positive' and 'negative' 
polynomials as before, now setting 7 =1= 6 = — n, yields precisely the equations (p.31|) and 
T32D. 



Regimes 3+ and 4+ 

As explained previously, the form of the solution of the recurrence relation for regimes 
and 4"*" can simply be obtained by replacing q with 1/g in ( |5.36D and ( |5.26| ), respectively. It 
should again be stressed that the precise meaning of the variable q in the various regimes is 
not related, and is given by equation ( |5.3| ) in regimes 1"^ and 2"^, and by ( 5.12| )in regimes 
and 4+. 



5.4 Thermodynamic limit 

In this section, using the solutions for the one- dimensional configuration sums X^'^ and 
y^'"^, we obtain expressions for the local height probabilities P^^(a). 



Regime 1+ 

As described in section E?T], in regime l"*" we have only ferromagnetic groundstate configura- 



tions. Hence we need to calculate 

m^ooj2^^^^q-a^x/ns^a)X^''\q) 



b = l,3,...,/,/ + l,/ + 3,...,L- 1. 



(5.39) 



From the solution ( |5.26| ) for the one-dimensional configuration sums, we find that we 
have to consider the m ^ 00 limit of the auxiliary function defined in ( ^.23|) . To do so, 
we need the result 



lim q 



k{k+a) 



k=—oo 



m 
k,k + a 



Qiq) 



(5.40) 



which holds for arbitrary fixed a. To establish this we take the limit inside the sum and use 
the elementary result 



lim 

m— >oo 



to obtain 



m 

k,k + a 



{q)k{q)k+a{q) 



1 



m—2k—a 



{q)k{q)k+a 
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k{k+a) 



fc=— 00 



m 
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00 k{k+a) 



q k=0 



iq)k iq)k+a Qiq) 



(5.41) 



(5.42) 



Here the last equality follows from the g-analogue of Rummer's theorem [|T9 



k=0 



(1 - g) ... (1 - g^)(l - z){l -zq)...{l- zq^-^) 



fc=0 



zq" 



(5.43) 



by setting z = g""*"^. From the above considerations we conclude that, in the m — oo hmit, 
the auxihary function yields the Virasoro characters defined in equation ( p..4|) 



|^(L+2)(L+l)j2 + [(L+2)a-(L+l)s]i_^(L+2)(L+l)i2 + [(L+2)a+(L+l)s]i+as| 



QiQ) j=-oo 



q2 

a'^X/n ^(i+2) 



(5.44) 



{<!)■ 



In the last step we have used A = ttL/{A{L + 1)), and have omitted terms independent of a. 

Substituting the appropriate elements of the solution ( ^.26|) , the X^^, into the expression 
for local height probabilities, and using the limiting behaviour of F^, we find 



where 



h 6 = 1,3,...,/ 

6+1 5=/+i/+3 



L-l. 



(5.45) 



(5.46) 



We note that s takes the values 1, 3, . . . , L. 



By performing the conjugate modulus transformation ( |C.8| ) we rewrite the Virasoro char- 
acters as 



-c/24-l/4h(/i-l) 



Q{q) 



71 



1/2 



:^3(f te-t),t)-^3(f fe + f),t)] (5.47) 



where 

As a result the local height probabilities can be written in the form 



Ni{s) 



(5.48) 



(5.49) 



where the normalisation factor in the denominator is given by 



a=l 
2L+1 



a=0 



L+2 J ' 



t). 



(5.50) 



We recall that the nome q is defined in terms of the nome p as q = exp(— 127rA/e). This 
yields the following relation between the nomes p and t: 



p 



^3L(L+2) 



(5.51) 



From the definition of tlie crossing factor 



and the simple identity 



(5.53) 



we find an alternative form for S: 

where we have again neglected a-independent factors. Substituting this back into the def- 
inition of the normalisation factor, we find that the summation in Ni{s) can actually be 
performed to give 



Nl{s) ~ 2{L + 1)Q (tl2L(L+l)) Q (^12L(L+2)) (^E^ ^2(L+l)(L+2) ) (^gzL^ t6L{L+l) 



(5.55) 



The proof of this denominator identity, being rather technical and lengthy, is given in ap- 
pendix [B.l| . 

Regime 2+ 

The working for regime 2"*" is very similar to that of regime l"*". Again we have only ferro- 
magnetic groundstates, and we are interested in calculating 

P'»=J^3. J .-.JJ.rJ'L.. 6=1,3,. ..,/-l,/ + 2,/ + 4,...,L-l. (5.56) 



From the solution (|5.36|) for the one-dimensional configuration sums we see that we have 



to take the limit m ^ oo of the auxiliary function C"^ defined in (|5.33| ). Using the result 
( CT) yields 

„4(a— r)a oo 

lim = 2_ V fn^L+l)Lf~[{L+l)r-La]j_ {L+l)Lf + [{L+l)r+La]j+ar\ 

"^^^ '"^ ' Qiq) ^ 

^ ^i(a~.)a-A(f;+^)+c/24 ^(^+1) (^) (5.57) 



where we have used that A = 7i{L + 2)/(4(L + 1)). 

Substituting the elements Y^^^ of the solution (|5.36|) , and using the limiting behaviour 



of Gin, ^6 find that the local height probabilities are given by 

Ea=lS{a)Xr,a ' [q) 



where 



b 6=1,3,...,/-! 

6-1 6 = / + 2,/ + 4, 1. 



(5.59) 



We note that r takes the values 1, 3, . . . , L — 2. 

If we use the conjugate modulus form of the Virasoro characters, we can rewrite this as 

S{a] 



N2{r) 



;^3(f (i-iTT)'0-^3(f (i + ZTi)'^ 



where the normalisation factor in the denominator is 



2L+1 



a=0 



(5.60) 



(5.61) 



From the relation between the nomes p and q we find that p and t are again related as in 
equation (|5.51|) . Inserting the form ( ^.54|) for the crossing factor S into the expression for 
the normalisation factor A'"2, we can carry out the summation over a. The result, proved in 



appendix p.2| , is 



f-{L+2){2L+l)/2 ^^^^^ ^6L(L+2)' 
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^^(^2r^^^6(L+l)(L+2)j 



(5.62) 



Regime 3"*" 

The solution of the recurrence relation for regime 3"'" is obtained by replacing g by 1/g in the 
solution ( p.36| ) for regime 2"*". We therefore need to consider the auxiliary function GJ^, with 
q replaced by 1/q. The effect of this replacement on the Gaussian multinomials is given by 



m 
kJ 



k'^+fi+kl-{k+l)m 



1/9 



m 
kJ 



(5.63) 



Applying this to G!^, we find 



ml 
oo 



X g 



^^-iL+l)Lf + [(L+l)r—La]j+k'^+[k+2{L+l)j+a-b]'^-m[2k+2{L+l)j+a-b] 



k,k+2{L+l)j+a-b 
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(5.64) 

}■ 



In the first term we replace by — A; — (L + l)j + (m — /i^ — a + 6)/2 and in the second term 
we replace k by —k — (L + l)j + {m — fia — a — b)/2, where = 0, 1 is determined by the 
requirement that 



(m - /ia - a + 6)/2 e Z. 



(5.65) 



The subscript of /i, indicating its dependence on a, is included for later convenience. After 
simplification, we thus obtain 



j,k=—oo 
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(5.66) 

}■ 



Antiferromagnetic phases In contrast to the previous two regimes, we now have an- 
tiferromagnetic as well as ferromagnetic phases. At first, we restrict our attention to the 
antiferromagnetic groundstates, and calculate 



P 



bb+l. 



q' 



1,2,4,6,...,L-1. 



(5.67) 



It follows from the solution ( |5.36| ) that we have to find expressions for the function GJ^ of 
equation ( p.66| ) in the m oo limit. To do so we use the result 
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+b 



(5.68) 



We therefore conclude that, in the limit m — oo, the auxiliary function ^ is given by a 
product of Virasoro characters 



lim q'^'/^G^^-'ib) = 
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-as)/2 - 



Qiq) t^o (q) 
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(5.69) 



x,Mq) xi^"\q) 



with A = 7r(L + 2)/(4(L + 1)). In the last step we have used the following Rogers- Ramanuj an 
identity [^] for the c = | characters Xo and Xi/2'- 



oo 2k^+2^iak 



fc=0 
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2k + IJ,a 



(5.70) 



In contrast to the ferromagnetic phases treated so far, the local height probabilities for 
antiferromagnetic phases depend on whether m is taken to infinity through odd or even 
values. We therefore write P^^"'"^(a), where a is defined to be the parity of m + 6 



a = m + b mod 2. 



(5.71) 



For fia this gives 

fia = cr + a mod 2. 



(5.72) 



Replacing g by 1/g in the solution ( |5.36| ) for the one-dimensional configuration sums and 
using the result (|5.69|) , we find 



bb+l. 



where 



h 

h + l 



6 = 2,4,...,L-1. 



(5.73) 



(5.74) 



We note that s takes the values 1, 3, . . . , L. 



Performing a conjugate modulus transformation, we find that the local height probabili- 
ties are given by 
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where the normalisation factor is defined as 
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(5.75) 



(5.76) 



From the relation between the nome p = exp(— e) and the nome q = exp(— 47r(7r — 3A)/e) we 
obtain 



p = t 



(5.77) 



We take the conjugate modulus transformation (|C.8|) of the c = | characters after first 
rewriting them using the formula PT| 



Xf.a/2iq) 



„Ma/2-l/48 



(5.78) 



Furthermore, we rewrite the crossing factor S by using the identity (|5.53| ) as well as the 
relation between the nomes p and t. Substituting this into the defining relation of the 
normalisation factor N^, we can perform the summation, yielding 



N,(s) 



2(L + 1)Q (^4(^-^-4)) Q (t2(L+l)(L+2)) 



X 



2{L+1)' 



^3 (f , (;^, t^L-2)(L+l 

The proof of the denominator identity is presented in appendix |B.3 



(5.79) 



Ferromagnetic phases For the ferromagnetic phases of regime the local height prob- 
abilities are given by 



P''Ha) = lim 



6 = 2, 4,..., /,/ + !,/ + 3,..., L. (5.80) 



From the expressions for Y^^^ in equation (|5.36|) it follows that we have to consider the 
following combinations of the function GJ„ in the m ^ oo limit: 



G'^\h)-q-^-'G'^\{b+l) 



6 = 2,4,...,/ 



6 = / + l,/ + 3,...,L. 



Here is given by equation ( |5.66|) and we have neglected terms that have an extra factor 
g™, which do not contribute in the large m limit. We again take m oo through even or 
odd values of m with a and h fixed and /i^ = 0, 1 chosen accordingly. For this we need the 
result 



lim q 



2m 

m — a — b,m — a + b 
= lim 



2m - 1 
m — a — b — l,m — a + b 



l-q 



2m 



2m 

m — a — b,m — a + b 



— lim q 



-a—b 



2m 

m — a — b,m — a + b 



(5.81) 



where we have used equation (|5.68|) . If we apply this limit to the combination of 
functions, it follows that 



^im^g('"'-'"+^)/2 \pb^\h) - q-^^^Gl+l.ib + 1) 

5,2 _ 



Q 



1 oo ^2fc2 + (2^„-l)fc 

(b-a)b/2 ^ y 



oo 

X > ^q 

j=~oo 



|^(L+2)(L+l)i2 + [(L+2)a-(L+l)fe]j _ ^(L+2)(L+l)i2 + [(L+2)a+(L+l)b]j+afe j (^^,82) 



= y^^-'^)^/^-^-'"^/"-^/"xi/i6(g)xir(g) 

To obtain this result, we have used the Rogers- Ramanuj an identity for the c = | character 



oo ^2k^ + {2na-l)k Q[q^) 



E 



Q{q) 



(5.83) 



The left-hand side of this identity does not depend on the value of /Iq, and hence the result 
( ^.82|) is independent of the parity of m. For a ferromagnetic phase this is indeed what one 
would expect. Similarly we find that 



lim g(™^-™-'')/2 



Gl{h)-q-'"^^'G\;Uh-l 



(5.84) 



Replacing g by 1/g in the solution ( |5.36| ) for Y^'"^, and substituting the above results, 
we get 



Y.'a=iS{a)xiT\q) 



where 



_{ h 6 = 2,4,...,/ 

^~\6 + l 6 = / + l,/ + 3,...,L. 

We note that s takes the values 2, 4, . . . , L + 1. 

If we continue as for the antiferromagnetic phases, we arrive at the final result 

S{a) 



:^3(f (iTT 



L+2 I ' 



t 



■^3(1 {-J^,+ 



L+2 



where the term occurring in the denominator reads 

2L+1 



a=0 



,L+1 L+2 J ' 

Again the summation over a can be carried out, yielding 

Sia) t-^(^+2)/2 l^n^ t(L+l)(L+2)\^ t{L-2){L+l) 



4(L + l)g2 (t4(L+l)(L+2)) Q {t^iL'^^i)) 

'^IV^+l' / \^ 2(L+1) ' ' 



anL 



-L^-i 



X 



di (;^, t2(L-2)(L+l)^ ^3 (^^^ t2(L-2){L+l 

A proof of the denominator identity leading to this result is given in appendix |B.3 



(5.85) 



(5.86) 



(5.87) 



(5.88) 



(5.89) 



Regime 4+ 

The solution of the recurrence relation for regime 4"*" is obtained by replacing g by 1/g in 
the solution ( [5. 261 ) for regime l"*". Consequently we need to consider the auxiliary function 
F^, with g replaced with 1/g. Using the inversion formula (|5.63| ) and carrying out the same 
sequence of transformations as for regime 3"*", we obtain 



:p2b-r 
m 



{h) = q 



{b'^+lia-m?-ar)/2 



^ ^2k^+2fj.ak 

j,k= — OD 



X i^q{L+l)Lj^-l(L+l)r-La]j J 



{-m-fia-a+b)/2-k-{L+l)j,(m-fj.a+a-b)/2-k+{L+l)j 



(5.90) 



_q{L+l)Lj^ + [{L+l)r+La]j+ar 



{m-tMa-a-b)/2-k-{L+l)j,{m-fj.a+a+b)/2-k+(L+l)j 



Antiferromagnetic phases We again have antiferromagnetic as weh as ferromagnetic 
phases. Both types admit treatment akin to that apphed in regime S"*". In the following 
we therefore leave out some of the details. We begin by treating the antiferromagnetic 
groundstates and calculate 



P^'+^a) = \im ' .^ 6 = 1,3..., L- 2. (5.91) 

From the solution (15.261) we see that we have to obtain the m — >■ oo limit of the function Fl, 



in equation (|5.9CI|) . We find that in this limit we get a product of Virasoro characters 

Jim g-V2^2.-.(^) ^ g-'^'^/-X,„/2(g) Xi^^'\q) 
where A = 7rL/(4(L + 1)). For the local height probabilities this gives 

pbb+l . ^ _ 'g(Q)XM./2(g)X^^a+'^(g) , . 

Ea=iS{a)Xt^^/2{q)Xr,a ' [q) 

where 

r = b 6 = 1,3,..., L- 2. (5.93) 
After a conjugate modulus transformation this gives 

:^3(f (i-j^) ,t)-^3(f (i + j^) ,t)] (5.94) 



with the following normalisation factor: 

2L+1 

N,ir)= Y: ^(«)x.„/2^3(f (i-xfr)'^)- (5-95) 

a=0 

From the relation between the nome p and q we get 

j9 = t^(^+^). (5.96) 
Using the result of appendix [B.4] for the normalisation factor A5, we finally obtain 

(5.97) 



j^-L(L+3)/4^ ^il+2fiah j^L(L+l)/4j 


^,(^^,t2L(L+4)^ 




2(L+l)g (t4i(i+4)) Q (t2L(L+l)) ^3(|^tL(L+l)j 


^9^(^I^,t(i+l)(L+4)^ 





Ferromagnetic phases For the ferromagnetic phases of regime 4"'", the local height prob- 
abilities are 

P'^a) = lim f ^^Iw "^^i"''^?..'^ & = 2,4,...,/-l,/ + 2,/ + 4,...,L. (5.98) 

From the solution for the configuration sums X^'', as listed in equation ( |5.26| ), it follows 
that we have to take the the m — 00 limit of the following combinations of the function F^: 

F'^\h) - q~^^'F'^\{h -I) 6 = 2, 4, 1 

Fl{h) - q-^-'F'^\{h + 1) 6 = / + 2, / + 4, . . . , L. (5.99) 



where is given by equation ( p.90|) . The infinite hmit can easily be taken using (|5.81|) and 
( ^.83|) , and we arrive at 



lim - q-^^^'Fll^zMb - 1)] ~ g-'^^/^xIf.^'^lg) 

hni - q-'--'F!^t\{b + 1)1 ~ g-'^'/^X&a 



(5.100) 



Again we observe that the dependence on the parity of m has dropped out. 

Replacing g by 1/g in the solution (|5.26|) for and substituting the above results, we 

get 



P'^ia) 



where 



_ \ b 6 = 2,4,...,/-l 

^ ~ \ b-l 6 = / + 2,/ + 4,...,L. 

We note that r takes the values 2, 4, . . . , L — 1. 



(5.101) 



(5.102) 



Taking a conjugate modulus transformation, we obtain the final result 



S{a) 



with the following normalisation factor: 

2L+1 
a=0 



:^3(f (i-iTT)'^)-^3(f 



(5.103) 



(5.104) 



Performing the sum over a, using the denominator identity for Nq in appendix p.4| , yields 



iV6(r) 



4(L + 1)Q2 (t4L(L+l)) Q (t4L(L+4)) 



X 



+l)(i+4) 



(5.105) 



6 Order Parameters and Critical Exponents 

Following Huse 0, we define generalised order parameters in terms of the local height 
probabilities obtained in the previous section: 



L sin 



a=l sin(f^) 
where = 0, 1, 2, . . . , L — 1 and 



P±ia) 



1 r 



P^''+^{a)±P^'{a 



(6.1) 



(6.2) 



For A; = these order parameters are trivial because -Rq + ~ ^ -^o - ~ 0- -^^^ 
romagnetic phases there is no distinction between the two sublattice probabihties, and so 

Rf,- = 0- 

We are not able in general to perform the sum in ( |6.1|) . To determine the associated 
critical exponents we therefore expand the local height probabilities for the four regimes in 
powers of t, and hence of p. In this process the critical values of the local height probabilities 
may also obtained, and are given by 

P^a) = sin2(ffr) (6-3) 

in all regimes, and are independent of the boundary spins b and c. 

Using the exponents obtained via the free energy, we extract the conformal weights 
corresponding to the critical exponents of the generalised order parameters. 



6.1 OfF-critical perturbation 



We first identify the scaling field corresponding to the elliptic nome p. From the scaling 
relations 

1 - a . 1 



2 — a 1 + 

and the exponents a and S obtained in section ^, we find that the conformal weight of the 
perturbing field is 



AiL+2) 
^2,1 

^1,2 



L+4 
4(L+1) 

L-2 



4(L+1) 
(4) , . (L+2) _ 1 



AC*) I A V^- 
^2,1 + ^2,1 

A (4) , A (L- 
^2,1 + ^1,2 



I L+4 
2 4(L+1) 



2 4(L+1) 



regime 1^ 
regime 2^ 
regime 3^ 
regime 



(6.5) 



We therefore conclude that the nome p corresponds to a perturbation in the direction of the 
spinless operator 0p, with 



LiL+2) 
'l,2 



'L-\,L^\ 
'L-\,L~\ 



,(4) , (L+2) 

'2,1 <Pl~i,l+i 
.(4) AL+l) 

Ha (Pl-i,l-i 



regime 1^ 
regime 2^ 
regime 3^ 
regime 4^. 



(6.6) 



Here denotes the operator of the unitary minimal model with central charge 
6 



6.2 Order parameters for regime 1+ 

Because this is a ferromagnetic regime, the only non-zero order parameters are 



L 



s,u(ltglF) 



Rt+ ^ -R? = E \r/ p^-M (6.8) 

a=i sm(^^_^J 

where A; = 0, 1, 2, . . . , L — 1. To evaluate the associated critical exponents, we expand 
equation ( 5.55|) using the representation (|C.1|) of the "i^-functions as a series of trigonometric 
functions. In the t ^ limit this yields 

N,is) 2(L + l)tsin(^) ^ ^""^ ^'-'^ 

Expanding also the T^s-functions in the expression (|5.49| ) for the local height probabilities, 
and using a simple trigonometric identity for the difference of cosines, gives 



9 gij2 ( _5iZi_ I CO 



(6.10) 



Setting n = 1 we obtain the aforementioned result for the critical values of the local height 
probabilities ( |6.2|) . 

If we substitute the above expansion into the expression for the order parameters, then 
use the orthogonality relation 



2 ^ 



Y Esm(f^) sin(fH^) = [Sn-n.,2kiL+l) - Sn+m,2kiL+l)] (6.11 



~^ a=l k=~oo 

and interchange the order of summation, we find 



-| oo 

fif = J2 t'^+^+^"-(^+^)i'~^sin( ^"+^+\"^^^+^)1^" ) +0(t3^(^+^)) . (6.12) 



m=— oo 



sm 

To leading order we thus find 

Hence, for > 1, the order parameters vanish at criticality with a power law behaviour 

Rl^ ~ p'/^'^ (6.14) 

where the critical exponents Sk are given by 
3L(L + 2 



(A; + l)^ 



A;= 1,2,...,L-1. (6.15) 



Here we have used the fact that p = t^iCi+z) ^Y^e corresponding conformal weights 

Afc we use the scahng relation 

A, = (6.16) 

Ok 

From this relation we get 

^-il^l^-^--' *^1.2.....L-l, (6.7) 

Thus we have obtained all 'diagonal' weights of the Kac table ( |1.2| ). 

6.3 Order parameters for regime 2+ 

To obtain the generalised order parameters for regime 2+ we closely follow the previous 
working. The non-zero order parameters are again given by ( |6.8|) , with k = 0, . . . , L — 1. 
Expanding the results (|5.6CI|) and ( |5.62| ) for the local height probabilities, we find 



sinf^j 



m=—oo 



For k = 1, . . . , L — 2, the leading term is obtained for m = 0. When k = L — 1, however, the 
amplitude of this term vanishes and we need the next-to-leading term m = — 1. Consequently 
we get 

Vi-itC^+D^-i A; = l,...,L-2 

sin(f) 

R't, ~ 2cos(^) (6.19) 
where we have used the fact that r is odd. From this we can read off the exponents 



To find the corresponding conformal weights A^ we again use the scaling relation (|6.16| ). In 



addition to the diagonal weights of the Kac table this yields the conformal weight A2 1 



(L+l). 



{k + lf~l 

AL{L + 1) 



Afc+i^fc+i k — l,2,...,L 2 



Al-1 - 4^(^ + 1) -^2,1 • (6.21) 



6.4 Order parameters for regime 3^ 
Antiferromagnetic phases 

Using the identities 

and the result (|5.79|) , it follows from equation (|5.75|) that 



P. 



bb+l. 



(6.22) 



(6.23) 



For the order parameters (|6 . 1|) this yields the relation 



(6.24) 
(6.25) 

Consequently we can restrict our attention to R'^'X^- 

As for the regimes l"*" and 2"*" we first expand the normalisation factor (|5.79|) . We note, 
however, that we now need more terms in the expansion: 
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^^(L+l)(L+2)/4 j^bb+1^ 
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l + 2t^'-^cos(f^) +0(t2(^-2)(^+i))] .(6.26) 



Ns{s)^ (L + l)t sin(- 

Expanding the ^^s-functions appearing in the local height probabilities ( 5.75|) we arrive at 
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^r^-isin(l^) sin(|ni 
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j anL 
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+0{t 
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(6.27) 



If we apply the orthogonality relation ( |6.11| ), after substituting this expression in the defini- 
tion of the order parameters, we obtain 
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L+2J 
+t 



[fc+l+2m(L+l)]s7r 
L+2 



.=—00 

[fc+l~L+2m(L+l)]2+L2-5 ^-^ ^ [k+l-L+2m{L+l)]s7v 



sm 



L+2 



(6.28) 



_^^[fc+l+L+2m(L+l)]2+L2-5 ^-^ ^ [fc+l+L+2m(L+l)]g7r ^ | ^ q ^^2{L-2){L+iy 

Except for k = L — 1, the leading term comes from the first term within the summation, 
with m = 0. Defining the exponents 6k^+ as in equation ( p.l4| ), we can read off 

L-2. (6.29) 



(fc + 1)2 - 1 



fc=l,2. 



When k = L — 1 the first and the third term, with m = and m = — 1, respectively, are of 
the same order. Using the fact that s is odd we find that R^l^l_^_ vanishes. Because of the 
symmetry ( |6.24| ) this is as it should be. 



Using the relation (|6.25|) we readily obtain the exponents 



^'--= (L + 2)(L+lV4"+(L-tp-l *=1.2.-.i-2. (6-30) 



From the scaling relation ( |6.16|) and the conformal weight Ap in equation (|6.5| ), we thus find 

A _ (k+lf-l _ A (4) , A (L+2) 

- 4{L+2){L+1) ~ ^1,1 ~^ ^k+l,k+l 



_ 1 , (L-kf-1 _ . (4) . (L+2) 
^k.- — iR "T ^/r I o\/r I 1^ — ^2.2 T 



A;= l,2,...,L-2. (6.31) 



^k- - 16 4{L+2)(L+1) ~ ^2,2 ^L-k,L-k 

We conclude that we have obtained all but one of the diagonal weights of the relevant Kac 
table, shifted by the diagonal c = | weights. It is this structure, which we observe also 
in Regime 4+ below, together with the products of Virasoro characters which we found in 
the local height probabilities for these two regimes, which led to the proposal ( p..7| ) for the 
modular invariant partition function for the corresponding critical branches. 

Ferromagnetic phases 

The derivation of the order parameters for the ferromagnetic phases of regime proceeds 
along similar lines to the calculation for the antiferromagnetic phases. 

The non-zero order parameters are again given by ( |6.8|) and, using the results (|5.87 ) and 
( p. 891 ), are found to be 



1 oo 

j^bb ^ ^ ( ,\k+l+2m(L+l)]'2-l ^-^ / [A:+l+2m(L+l)]s7r 



m=— oo 



^[fc+l-L+2m(L+l)1^+L^-5 ^ [fc+l-L+2m(L+l)]s7r j 



_j_^[A:+l+L+2m(L+l)]2+L2-5 gjn ^ [k+l+L+2m{L+l)]s7T ^ | _^ ^ ^^2{L-2){L+1) ^ 

For k = 1, . . . , L — 2 the leading term comes from the first term within the curly braces, with 
m = 0: 

V^±2; ^(.+1)2-1 fc = l,...,L-2. (6.33) 

For k = L — 1 the first term with m = and the third term with m = — 1 are of equal order. 
Using the fact that s is even gives 

i?ti~-4cos(^) t^'-\ (6.34) 

From these two equations we extract the exponents 

+ k = l,2,...,L-l. (6.35) 

which are the exponents 5k,+ found in ( |6.29[ ), plus the 'missing' exponent 



6.5 Order parameters for regime 4+ 
Antiferromagnetic phases 

From the resuhs (|5.75|) and (|5.971 ) for the local height probabilities, and the symmetry 
relation 



( H fL{L+l) 
^bb+l ^ '^H4''^ y ^66+1 



~ ^ti(i+i)/4^66+i (6.36) 
we obtain the following results for the order parameters in the small field limit: 

■ / (fc+l)r^ \ 



sm 



■ / (L-fc)r^ \ 

~ v^t^(^+i)/^ ^ /. ^ (6.37) 
sin(f) 

From this we read off the critical exponents 
L{L + A) 



5k, 



L{L + A) 

- LiL + l)/A + iL-ky-l ^^-^^^ 
where = 1, 2, . . . , L — 2. We apply the scaling relation ( |6.16| ) to the above results to find 

A, _ _ ^(4) ^(L+l) 

- 4{L+1)L ~ ^1,1 + ^fc+l,fc+l 



_ 1 (L-fc)2-l _ , (4) . (L+1) 
^fc - - 16 + 4(L+1)L - ^2,2 + ^L-k,L-k ) 



= l,2,...,L-2. (6.39) 



We have thus obtained the diagonal weights of the Kac table, with h = L + 1, shifted by the 
diagonal c = | weights. 

Ferromagnetic phases 

From the expression for the order parameters which survive in the ferromagnetic phases 
( |6.8|) , and the results for the local height probabilities ( 5.103|) and ( ^.105 ), we have 

g.^/ (fc+l)r^ \ 

Rk /!^\ t'^k+iy-' k = l,2,...,L-2. (6.40) 

sin I 



From this we obtain a subset of the exponents found for the antiferromagnetic phases. We 
are also able to determine that 



i?5f_^^_4cos(^j t(^+2r-i (6.41) 
which corresponds to the conformal weight A^^l + A^^^^\ 



7 The phase diagram 



The study of the solvable manifolds spanned by the elliptic nome and the spectral param- 
eter reveals some aspects of the role these manifolds play in a larger parameter space. In 
particular the off-critical regimes support coexistence between a specified number of phases, 
each associated with a groundstate, detailed in section [5.1|. The critical branches are the 



critical points terminating these coexistence lines. In the following we discuss some aspects 
of the phase diagram of the dilute A models. The discussion is not limited to the solvable 
manifolds, but is restricted to the region where the weights are ^2 symmetric. 

When the ^2 symmetry + 1 — ais obeyed, the weights permit a duality trans- 

formation. This transformation is a direct generalisation of the orbifold duality, relating the 
A2k+i and Dk+2 Temperley-Lieb models [Q, to the dilute A and D models. Since the A3 
and D3 Dynkin diagrams are identical, we have a dual symmetry of the phase diagram in 
the case L = 3. Because the phase diagram for L > 3 is complicated by the large number of 
states of the model, and by the absence of a dual symmetry, we limit detailed discussion to 
the case L = 3. Only afterwards do we mention which conclusions are valid for general L. 

The duality transformation for the case L = 3 can be performed by separating each 
degree of freedom into an Ising-like variable a = ±1, that discriminates between the states 
a = 1 and a = 3, and the variable s = 0, 1, that decides whether a state is a = 2, so 
that a = 2 + sa. For a given configuration of s, the Ising variables live on the sites where 
s = 1. Since these Ising variables take the same value on neighbouring sites, there is 
effectively only one Ising variable in each cluster of s = 1 sites that are mutually connected 
(directly or indirectly) by nearest neighbour links. Such clusters form an irregular lattice, 
with interactions between them wherever they are linked via a second neighbour bond. On 
this lattice an ordinary Kramers- Wannier duality transformation can be performed which 
results in new Ising variables on the dual lattice, formed in the same way, but by the clusters 
where s = 0. The s-variables, fixed in the procedure so far, are now each replaced by 1 — s. 
The resulting duality transformation on the weights is 





(7.1) 




where the weights are symmetric under exchange of the states 1 and 3, and the transformation 
rules are invariant under reflection in the diagonals of the elementary squares. 

The self-dual subspace contains all four solvable critical branches. It constitutes a tran- 
sition between a phase symmetric under exchange of the states 1 and 3 and one in which 
this symmetry is spontaneously broken. The nature of this phase transition is not the same 
everywhere in the self-dual subspace. The phase structure within this subspace, which is 
equivalent to the phase diagram of the 0{n) loop model is indicated in Fig. ^ It may 
be parametrised in terms of the weights pi, . . . , pg. If we keep the weights p2, • • • , fixed, 
we may describe the self-dual parameter space in terms of the weight pi of the ferromagnetic 
states and the weights ps and pg of the antiferromagnetic states. 

When pi as well as ps and pg are small, the transition (across the self-dual space) is of 
second order, in the Ising universality class. It is described by the solvable critical point of 
branch 2. The character of the phases is ferromagnetic. When the weight pi is sufficiently 
large the transition is of first order. These regions are separated by an Ising tricritical 
transition governed by branch 1. 

When, on the other hand, ps and pg are sufficiently large and pi is relatively small, 
antiferromagnetic configurations, in which one of the sublattices takes the value a = 2 while 
the other sublattice fiuctuates between a = 1 and a = 3, are favoured. There will be a 
phase separation between regions where one or the other sublattice has a = 2. The onset 
of this phase separation is an ordinary Ising transition, which appears to be independent of 
the transition across the self-dual subspace. It is described by the critical point of branch 
4. The independence of these transitions is refiected in the product structure of the order 
parameters on this branch, see e.g., equation ( p.92| ), and in the value of the central charge 



c = 1, see Table |. It implies that the transition across the self-dual subspace is of the same 
(Ising) universality class on either side of the Ising transition within the subspace. 

When the ferromagnetic weight pi and antiferromagnetic weights ps, pg are all large 
relative to the other weights, the system must make a choice between the first-order regime, in 
which the ferromagnetic configurations dominate, and the critical antiferromagnetic regime. 
From the fact that the surface tension between these phases can be made arbitrarily large, 
we conclude that the transition between the two regimes is first order. The point where this 
first order transition is met by the other two transitions within the subspace is the critical 
point of branch 3. 

We further observe that the transition between the ferromagnetic and antiferromagnetic 
region of the self-dual subspace, persists away from the self-dual subspace, see Fig. In 
fact, when the weights are taken far from self-duality, the model may be reduced to the 
interacting hard square model |^ ^ |2^, |2^, in which the particles are represented by 



a = 1 or a = 3, and the empty sites by a = 2. Of course, by dual symmetry, the roles are 
reversed in the opposite extreme. The hard-square model is known to have an Ising critical 
and tricritical transition into a sublattice-ordered phase. We expect the critical transition 
to join continuously with branch 4 and the tricritical transition with branch 3. 

We now see that the critical point of branch 3 plays a central role. It sits at the intersec- 
tion of two phase-transition manifolds, i.e., the ferro-antiferromagnetic transition and the 
self-dual subspace. Furthermore, it sits at the line where both transition sheets turn first 



order. 

For general L we expect the same topology of the phase diagram. Though there is no 
longer a dual symmetry that maps the weight space into itself, the subspace given by ( |2.2| ) 
still plays the role of a phase transition. The universality class of the transition across 
this subspace varies with L, and is that of regime III-IV of the L-state ABF model for 
branch 2, and of the L + l-state ABF model for branch 1. The transition between ferro- 
and antiferromagnetic regions remains Ising-like for all L. Where the transition manifolds 
intersect, the universality class is that of a direct product of an Ising and an ABF model 
(except where they are first order). Again, for general L there are directions in the parameter 
space in which the model can be reduced to the interacting hard square model. 

To conclude we remark that with the phase diagram proposed above we do not do full 
justice to the fact that the critical behaviour of branch 3 is that of the product of an Ising 
model and branch 1. 
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Figure 3: The phase diagram for L = 3. (a) The self-dual subspace. The encircled numbers 
represent the four solvable critical branches; branches 1 and 4 describe the transition of the 
entire curves on which they are located, and branch 2 describes the transition across the 
whole lower left section of the diagram. The double curve marks a first-order transition, (b) 
The phase diagram with an additional parameter that reverses sign under duality. In the 
extreme limits (front and back) the model reduces to the interacting hard square model. 



8 An Ising model in a field 

Because the model in regime 2"*" is in the universality class of the Ising model in a magnetic 
field, we take a special interest in this case, and collect the results in this section. 

As we have already observed, we can make the identification {1, 2, 3} = {+, 0, — } to relate 
the heights of the dilute A3 model to the states of a spin-1 Ising model. This particular model 
has the restriction that a '+' and a ' — ' spin may not be adjacent on the lattice. Because the 
nome p breaks the up-down symmetry, it plays the role of the magnetic field. In Regime 2"*" 
we see from Table | that the central charge of the model is c = |. 

Using result ([4.3|) we can rewrite the free energy as 



o /""^ coshQvra; coshSvrx smhlQuux sinh(157r — 16u)a; ■i2iekx 1 

inKitt) = z y I : — - — —— — — e dx. lo./i 

' X smh loTTX cosh lovra; 



The magnetisation of this Ising model can be expressed as a hnear combination of the 
local state probabilities: 



m = {+)- (-) = P^\l) - P'\3) = RY/V2. (8.3) 

Analogously, we can define a density for the spin-1 model, which is the probability that a 
site is occupied by a spin of either sign, 

p^{+) + (-) = P'\l) + P'\3) = i(l + Rl'). (8.4) 
When we group the following combination of elliptic functions: 

^9i(f,t4) Ar2(r = 1) 4v^ ' ' 
and make some simplifications, the magnetisation and density are 

m = 4V3t'K {E t''^) E (-t^O^ ^240j _ ^24^ |^_^18^ ^144^ ^ (^_^30^ ^240^ | 

p = 2V6K{^s{0,f')^s{0,t''')-d2{0,t'')^2{0,t''')}. (8.6) 
The leading order behaviour is 



In Ksing ~ 


pl6/15 






= pl/15 


(8.7) 


p-1 ~ 


= p8/15 





The first two expressions both lead to the exponent S = 15. The last exponent is related by 
scaling to the conformal weights of the energy and spin operators: Ai^2/(1 — ^2,1) = 8/15. 

9 Summary and discussion 

In this paper we have considered the infinite hierarchy of dilute A models. These models, 
which belong to the family of dilute A-D-E models, can be viewed as spin-1 generalisations 
of the solid-on-solid models of Andrews, Baxter and Forrester For any integer L > 2, the 
dilute A model defines a L-state solid-on-solid model, with heights labelled by the Dynkin 
diagram of the classical Lie algebra A^. At the critical point, the Boltzmann weights p.l|) 



of this solid-on-solid model obey the Z2 symmetry of the underlying Dynkin diagram, but 
away from the critical point this symmetry is broken for odd values of L. 

For each value of L the model has four different critical branches. Two of these branches 
provide new realisations of the unitary minimal series and the two other branches can be 
viewed as the direct product of this same series and a c = ^ model. Away from criticality 
the four branches yield eight distinct regimes. 

For all regimes we have calculated the free energy. From this we have obtained critical 
exponents a, for even L, and exponents 6, for odd L. Among these results is the magnetic 
exponent 5 = 15 of the Ising model. 



The main part of this paper is concerned, however, with the calculation of the gener- 
alised order parameters for the symmetry-breaking models of the dilute A hierarchy. This 
calculation involved the evaluation of new sums-of-pro ducts identities for theta functions of 
the Rogers- Ramanuj an type. From the resulting expressions for the order parameters we 
extracted a set of critical exponents 5k- For the models corresponding to the unitary min- 
imal series, these exponents correspond to the diagonal conformal weights of the relevant 
Kac table. For the models corresponding to a direct product of a unitary minimal model 
and a c = I model, the critical exponents yield the diagonal weights of the Kac table of the 
minimal model incremented by either one of the diagonal weights of the c = | model. 

We have not yet succeeded in computing the order parameters for the non-symmetry- 
breaking dilute A models, obtained for even values of L. In this case the Boltzmann weights 
no longer enjoy the diagonal property ( |5.5|) and an additional diagonalisation is required. 
We hope to report on these 'even' members of the dilute A hierarchy in a future publication. 

Another open problem is the relation between the symmetry-breaking spin-1 Ising model, 
corresponding to the dilute A3 models in the regimes 2"^ and 2~, and the integrable field 
theory of the critical Ising model in a magnetic field, found by Zamolodchikov p5|, . As 



shown by Smirnov |30[, Zamolodchikov's S-matrix corresponds to an RSOS projection of the 
Izergin-Korepin R-matrix. Hence the dilute A3 model provides a likely candidate for describ- 
ing the corresponding solvable lattice model. The precise relation between Zamolodchikov's 
^'-matrix, which has a structure related to the Lie algebra Eg, and the dilute A3 model 
remains, however, unclear. By studying the excitation spectrum of the dilute A models, and 
in particular the A3 case, we hope to establish the connection, if present, in the near future. 
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A The function H 

In this appendix we give the weights in a form which is more suitable for considering 
the ordered or strong- field limit, and hence for obtaining the weight function H defined in 
equation ( p.5|) . 

We first define the new variables 



w 



^-2nu/e 2,_g-2^AA_ ^^^^^ 



In terms of these variables, the ordered limit is given by a; ^ with w fixed. Using the 
conjugate nome expressions ( |(J.^ ) of the theta functions and setting 

e{u,x^/^) = E{u) (A.2) 

the Boltzmann weights in conjugate modulus parametrisation read 



a a 



E(x^w~'^)E(x^w) ^E(w)E(x^w-^ 



a a / Eix^)Eix^) E(x^)E(x^] 



^(a±l a\^fa a \ 9a ^E{x^w-')E{-x^^-+'w-') 



a a J ya a±l J g^^^ E{x'^)E{-x^^<^ 

w( I, ^Ui^f^ 

\a±laj \a a j 



-w ^x" I —x^^ 



Qa ^ - I - ^(_a;2{a±l))^(^3,4a) ) E{x^)E{-X^^''+^ 



\aa±l/ \a a 



^[ ^2(_3,±2a+l) ) E{X^)E{X^) ^ 



yy(a±l a \ _ 9l ^ E{x^w-^)E{x^w-^) 
V « ga+iga-i"^ E{X^)E{X^) 



a 



a 



a±l a J gl E(x3)E(x±4a+2j 



gl E(-x2("±i))E(x4«)E(a;3)E(a;±4«+2) 
fail ^-1 ^(a;^^)^(x±^''-%) gl^^ E{w)E{x^^''-^w) 



X x^^'"—^, J , J + X' 



where we have omitted an overall normalisation factor exp(2-u(3A — M)/e) and ga is defined in 
equation ( ^.61 ). In this alternative representation of the Boltzmann weights, we can readily 
extract their leading behaviour in the ordered limit. To do so, we make use of the simple 
properties of the function E as listed in equation (|C.4| ) of appendix p. It is clear from the 
definition of x that the result depends on the value of A. For regime 1^, using definition 
( ^■251) , we get 
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(A.4) 



ff(6- 1,6,6- 1) 
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1 6 G {2,3, 1} . 



To facilitate the proof of solution ( p. 261 ) we include the following two terms to the above list: 
H{0, 1, 1) = and H{L + 1, L, L) = i(L + 1). 
For regime 2"'", using definition ( |5.35|) , we get 
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(A.5) 



iy(6+l,6,6+l) 
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r -ib-i) 
[ -ip-l) 

H{b±l,b,bTl) = 1 & G {2,...,L- 1} . 
Now we only add the single term H[0, 1, 1) = 0. 

B Denominator identities 

In this appendix we prove the six denominator identities used in section to normalise 
the local height probabilities. All identities involve sums over products of elliptic functions 
similar to the sums-of-products identities of Rogers and Ramanujan. 

For brevity we shall adopt the convention that sums over j, k,l,m and n always run over 

Z. 



beh\{l} 
b e h 
beh 
b eti 



B.l Regime 1+ 

For regime l"*" the denominator identity reads 

2L+1 



a=0 

( 

L+2' 



2{L + 1) t^(2L+3)/2g |^^12L(L+2)^ g (il2L(L+l) j (^.^ i2(L+l)(L+2)^ 



where s = 1, 3, 5, . . . , L and L = 3, 5, 7, . . .. 

Proof: To prove this identity we start by recasting the the left-hand side using the 
representation ( |(J.1D of the t9-functions as infinite sums: 

2L+1 

= _ i ^3L(L+2)/2 Y^_^y ^^ia{2j+l+kL+e)/{L+l) ^^7Tise/{L+2) 

a=0 j,k,i 

^ ^6L(L+2)(j2+j)+3L(L+2)fc2+<?2^ ^^ 2) 

Summing the roots of unity over a, we find that the above expression vanishes unless i = 
—2j — 1 — kL + 2m(L + 1), with m G Z. Hence we get 

Ni = -2i{L + 1) t3^'(i+2)/2 ^_'^y+k ^7Tisi2j+l-2k+2m)/{L+2) 

j,k,Tn 



X 



^6L(L+2)(j2+j)+3l,{L+2)fe2 + [2j+l+fcL-2m(L+l)]2 g-j 



where we have used the fact that s is odd. 

We now carry out a sequence of transformations to diagonalise the quadratic exponent 
of t. First we make the substitution m ^ m — j + k followed by — > A; + 2j. Then we 
replace the sum over A; by a sum over i and a by setting k = 3i + a, where a has to be 



summed over 0, 1 and 2. Similarly, we replace the sum over m by a sum over n and (3 by 
setting m = 3n + (3. Finally we make the replacements j j — i and £ i — n. After all 
these transformations we obtain 
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Ni = -2i{L + 1) t(3^'+6-f^+2)/2 J2 (-1)^'+" Q^is(6n+2l3+l)/{L+2) 

Q,/3=0 j,i,n 

.4(a+f3){a+l3-l)+2L(5a^+4f3^+6al3-a-2l3)+4L'^(a^+t3^+al3) 



^ ^6L(L+2)(3j+l+2a)j+6(L+l)(L+2)(3^-l+2«+2/3)^+6L(L+l){3n+l+2/3)n 

-2i(L + l)t(3i^+6^+2)/2 ^B_4) 

2 

^4(Q+/3)(a+/3-l)+2L(5a2+4/32+6a/3-Q-2/3)+4L2(Q2+^2_,_Q,^-) 



2 

X 

a,/3=0 



^ g7ris(2/3+l)/(L+2)_g^_g67ris/(L+2) ^12L(L+l){2+/3) ^ ^36L(L+1)^ 

X £'|^il2L(L+2)(2+a) ^36L(L+2)^ ^ |^_^12(L+l)(L+2)(l+a+/3) ^36(L+l)(L+2)^ 

where, in the second step, we have used the definition ( |C.3| ) of the function E. The terms 
with a = 1 are identically zero and the two terms with f3 = 1, a ^ 1 cancel. The remaining 
four terms factorise to give 

iVi = 2 i (L + 1) t(3i'^+6i+2)/2 ^-nis/iL+2)Q |^^12L{L+2)^ 

X "^|^_^12(i+l)(i+2)^^36(L+l)(L+2)^ _ ^4{L+l)(L+2) ^ -f-36{L+l){L+2)J 
X "^(_e6--/{^+2) ^12L{L+1)^ ^36L(L+1)^ (3 5) 

_g27ris/{i+2)^|^_g-67ris/(L+2) ^12L{L+1) ^36L(L+l)y 

We now note the following two elliptic function identities: 

E {-x'p,p') - xE (-x-'py) = ^^Pjfy'f 

hi \ x,pj 

E [-p^,p^) -pE = e'^'/^E (e'^'/^p) . (B.6) 

To prove the first identity we note that the function /(x), defined as the ratio of the left-hand 
side over the right-hand side, has the following periodicity property: fipx) = f{x). Since 
the function E{x,p) is analytic in < |x| < oo, the only possible poles in the period annulus 
P < |x| < 1 are the zeros x = ±1 and x = ip^/^ of the function E{x'^,p). Using formula 
( |U.4| ) to manipulate the function E it is readily verified that these poles have zero residue. 
Hence, by Liouville's theorem, / is constant. Setting x = —p^/'^ shows that this constant is 



one. The second identity is actually a corollary of the first. Replacing p by p^, then setting 
X = and again using formula (|C.4|) , yields 



E - pE = -p ' > I 



E {—p'^,p^) 



n=l 



= e^'/3E (e^'/^p 

Applying both identities and the relation ( |(J.(j| ) of appendix |C|, we obtain the desired right- 
hand side of (PI ). 

B.2 Regime 2+ 

For regime 2"^ we have to prove the following identity: 



a=0 

= 2(L + l)t(^+2)(2L4-l)/2g Al2L(L+2)\ g A 12(L+l)(L+2) \ ^ /^^2L(L+1)\ ^l^^J ^ 



where r = 1, 3, 5, . . . , L — 2 and L = 3, 5, 7, . . .. 

Proof: The proof of this denominator identity is similar to that presented for regime l"*". 
Using the representation of the "i^-functions as infinite sums yields 

2L+1 

= -it^L(L+2)/2 ^(^_^yg7ria(2i+l+fcL-<?)/{L+l)g7rirf/L^6L(L+2)(j2+j)+3L{L+2)fc2+<?2 ^^g^g^| 

a=0 j,k,l 

Performing the sum over a gives £ = 2j + 1 + fcL + 2m (L + 1), where m assumes integer 
values. Thus we find 

j,k,m, 



X 



^6L(L+2)L(j2+j)+3L(L+2)fe2 + [2j+l+fcL+2m(L+l)]2 -j^q^j 



where we have used the fact that r is odd. 

As before, we carry out a sequence of transformations to diagonalise the quadratic ex- 
ponent of t. First we make the replacements m ^ m — j and k k + 2j. Then we set 
k = 3i + a and m = 3n + [3 and finally we substitute j ^ j — i followed hj i ^ i — n. After 
making these transformations we obtain 

2 

]\[2 = —2i{L + 1) ^{3^^+6^+2)72 Y (^_X)"t4'3('3+l)+2i{3a^+4/32+2a/3+a+2/3)+4L2{a2+/32+a/3) 

a,l3=0 

^ g7rir(2/3+l)/L^|^_g67rir/L^12{L+l){L+2){2+/3)^^36{L+l){L+2)^ (B.ll) 
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The terms with a = 1 vanish and the terms with (3 = 1 cancel. The remaining four terms 
can be written as 



fB.12) 



iV2 = 2 i (L + 1) t(3i^+6^+2)/2 e-^ir/Lg |^^12L(L+2) 

X ■^|^_^12L(L+1)^^36L(L+1)^ _ ^4L(L+1)^^_-|^^ ^36L(L+1) 
X "_/^^_g67rir/L^12{L+l){L+2)^^36{L+l){L+2)^ 

_g27rir/L^/ g-67rir/L ^12(L+l)(L+2) ^36(L+l)(L+2) 



If we again apply the two identities of equation ( [B.(j| ) and use the relation ( |U.(j| ), we obtain 
the required right-hand side of equation (p.8| ). 

B.3 Regime 3+ 
Antiferromagnetic phase 

For the antiferromagnetic phases of regime the denominator identity is given by 

2L+1 
a=0 

= 2(L + 1) t(^-l)(^+2)/4g |^^4{L2-4)j Q (^^2{L+l){L+2)^ (g^^^g) 

X^3(|,t(^+l)(^+2)^^^|^_s^^^(L+l)(L-2)^) 

where = a + o" mod 2 with a = 0, 1; s = 1, 3, 5, . . . , L and L = 3, 5, 7, . . .. 

Proof: We split the summation on the left-hand side depending on the parity of a. Using 
the fact that i)3{u + |vr,p) = {}4^{u,p) = ^3{u, — p), we get 
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We now represent the theta-functions with a-dependent arguments by series to obtain 
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Performing the sum over a, we find that i = —2j — 1 — kL + m{L + 1) and hence 
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where we have used the fact that s is odd. 

Again we carry out a sequence of transformations to diagonahse the quadratic exponent 
of t. First we make the replacements m —>■ m — 2j + 2k and k —>■ k + 2j, followed by the 
substitutions k = 3i + a and m = 4n + /3. Finally we set j —>■ j — i and i ^ i — n to obtain 

2 3 
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In this expression the terms with a = 1 vanish and the terms with /? = 1 and (3 = 3 cancel. 
Using the simple identity 

^3(f,p)+^3(f,p)=2i?3(f,/) 

which follows directly from (|C.1| ), the remaining four terms factorise as 
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We now note the identities 
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The first one follows directly from Liouville's theorem and the second one can be obtained 
from the first upon setting x = p after replacing p hj p^. If we apply these two identities, 
as well as equation ( |C.6| ), we find precisely the right-hand side of equation ( |B.13|) . This 
completes the proof. 



Ferromagnetic phase 

The denominator identity for the ferromagnetic phases in regime 3"*" reads 

2L+1 
a=0 

where s = 2, 4, 6, . . . , L + 1 and L = 3, 5, 7, 

Proof: We begin by writing the sum over the ■j?-functions in the left-hand side as 

2L+1 

N4 = _ ii(^^-4)/2 ^ E e'''"^^^'+^+''^+^^/^^+^^ g-7ris^/(L+2) ^(L2-4)(2j2+2j+fc2)+^2^ (B.22) 

a=0 j,k/ 

Summing over a yields i — —2j — 1 — kL + 2m{L + 1) and hence, using the fact that s is 
even, we get 

N4 = -2i(L + E (-l)^ e'''^(^^'+^"^''+^"')/(^+^) 

^ ^(L2_4)(2j2+2j+fe2)+[2j+l+feL-2m(L+l)]2_ (B.23) 

In order to diagonalise the quadratic exponent of t, we make the replacements m — > 
m — j + k and k ^ k + 2j followed by the substitutions k — S£ + a and m — 2n + p. Further 
replacing j ^ j — i and i ^ £ — n leads to the result 
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The terms with a = 1 vanish and the other terms may be combined to give 
7V4 = 2 i (L + 1) e--/(^+2) i(^'-2)/2Q (^4(^^-4)) 
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If we multiply this expression with the prefactors outside the sum in equation ( [B.21|) and 
use the following identities 

E ix p,p — xE [x p,p ] = V2e ' — 

\ / \ / E[i,p) 

E +pE [py^] = V2e-/^ 1^ (B.26) 

we find precisely the right-hand side of equation ( |B.21| ). The proof of these two relations is 
similar to those obtained earlier. The first identities follows from application of Liouville's 
theorem and the second identity is a consequence of the first. Replacing the nome p by p^ 
and then setting x = p^ yields 

E{p\p^-) +pE{p^p^') = V2e--/^g(/) ^^^^ 
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B.4 Regime 4+ 
Antiferromagnetic phase 

The denominator identity in the antiferromagnetic phase of regime 4"*" is given by 

2L+1 
a=0 

= 2(L + l)t^(^+3)/4g (^^4L(L+4)^) g (^^2L(L+1)^) ^ ^L(L+1)^ ^ t(i+l)(i+4)) (B.28) 

where /Iq = a + a mod 2, with a = 0, 1; r = 1, 3, 5, . . . , L — 2 and L = 3, 5, 7, . . .. 

Proof: Following the proof for the antiferromagnetic phase of regime 3"'", we split the 
summation on the left-hand side depending on the parity of a. Representing the theta 
functions with a-dependent arguments by series yields 

L 

jY^ = _ •^L{L+4)/2 ^27ria(2i+l+fcL-£)/{L+l) ^tt i rf/L ^L(L+4){2i2+2i+fc2)+£2 

a=0 j,k,e 



X 



Summing the roots of unity over a we see that i = 2j + 1 + kL + m{L + 1) and as a result 
we have 

iVg = - i(L + l)t^(^+'^)/2 ^ (^_]^y+feg7rir(2i+l+m)/L^L(L+4){2i2+2i+fc2)+[2i+l+fcL+m(L+l)]2 

X [(-1)^^93(1, (_l)-t^(i+l)/4j + ^3(|E^ (_l)-tMi+l)/4^] (B 30) 



where we have used the fact that r is odd. 

We diagonahse the quadratic exponent of t by carrying out the following sequence of 
replacements. We first set m — > m — 2j and k ^ k + 2j . We then replace k = 3i + a and 
m = in + P and finally we substitute j ^ j — i followed by £ — > £ — ra. We then have 

2 3 
a=0 /3=0 

X "(-l)/^^3(|,(-l)'^t^(^+i)/4) +^3(|:,(-l)-^^(^+l)/4y 
^ g7rir(/3+l)/L^|^_g47rir/L^2(L+l)(L+4)(3+/3)^^8(L+l)(L+4)j (B.31) 



X 



^|^^4L(L+4)(2+a) ^12L(L+4)^ ^ |^_^2L(L+l){5+4a+3/3) ^24L(L+1 



As before, the terms with a = 1 vanish and the terms with (3=1 and (3 = 3 cancel. Using 
the relation ( p.l8|) , the remaining four terms can be written as 



N, = 2i{L + l)e--''^/^t(^'+^^+2)/2g (^4L(L+4)^) ^3(^.^^L(L+1) 

y- ' ]^ ^_f-WL{L+l) ^24L(L+1)^ _ ^2L(L+1)^|^_^2L(L+1) ^24L(L+1) 

^ |^_g47r i r/L ^2(L+l)(L+4) ^ ^8(L+l)(L+4) ^ (B.32) 
_g27rir/L^|^_g-47rir/L ^2(L+l)(L+4) ^8{L+l){L+4) ' 



After application of the identities ( p.2CI|) , we obtain the required right-hand side of ( |B.28|) . 



Ferromagnetic phase 

The denominator identity for the ferromagnetic phase in regime 4+ is 



E ^i(ifT>^'"^"^'0^3(^,t^(^^^V3(f (l-ITl)'^) (B-33) 

a=0 

= 4(L + 1) t^(^+2)/2g2 (^4L(L+1)^ g (^4L(L+4)> 



^9^(l^^t2(L+l)(L+4)^ 


^3(!^,t2(L+l)(L+4)^ 




^9l(f,t(i+l)(^+4)) 





where L = 3, 5, 7, . . ., and r = 2, 4, 6, . . . , L — 1. 

Proof: We rewrite the sum over the '(9-functions in the left-hand side as 

2L+1 

jVg = _ i ^i(i+4)/2 ^ g7ria(2j+l+fcL-£)/(L+l)g7rirf/L^L{L+4)(2j2+2i+fc2)+£2_ (B.34) 

a=0 j,k,l 

By summing over the roots of unity, this vanishes unless £ = 2j + 1 + kL + 2m{L + 1). As a 
result we find 

Nq = -2i{L + (^_]_-)ig7rir(2j+l+2m)/L^L(L+4)(2j2+2j+fc2)+[2i+l+fcL+2m(L+l)]2 

j,k,m 



where we have used the fact that r is even. 

This expression is diagonahsed by carrying out the sequence of transformations m 
m — j,k^k + 2j,k = 3i + a, m = 2n + f3,j^j— i and i ^ i — n. This gives 

2 1 

= —2\(L + 1) ^(-^^^+4^+2)72 ^ ^ ^4/3(/3+l)+2L(2a2+4/32^2a/3+Q:+2/3)+2L2(a2+2/32+2a/3) 

a=0 p=0 



X 



X 



gTT i r(2/3+l)/L _g ^^471 i r/L ^2(L+l)(L+4)(3+2/3) ^ 
^^^4L(L+4)(2+a) ^12L(L+4)^ ^ |^^2L(L+l){5+4a+6/3) ^24L(L+1)^ 



(B.36) 



Here the terms with a = 1 vanish. The other terms may be combined to 
Ne = 2i(L + l)e-'^''^/^t(^'+^^+2)/2g |^^4L(l+4)^ 



^^^10L(L+l) ^24L(L+1)^ _^ ^2L(L+1)^|^^2L(L+1) ^24L(L+1) 
^^^inir/L ^2{L+l){L+4)^ ^8(L+l)(L+4)^ 

_Q2nir/L^^^-47Tir/L ^2(L+l)(L+4) ^ ^8(L+l)(L+4)^ 



(B.37) 



If we again use the identities (p.26| ), and multiply the above expression by the correct 
prefactors from the left-hand side of ( p.33|) , we indeed obtain the right-hand side of (p.33| ). 



C Theta functions 

In this appendix we list several definitions and relations for the Jacobian -(^-functions, used 



in the main text. For a more complete introduction, we refer the reader to e.g., pT 



The four standard -(^-functions, of nome p, \p\ < 1, are defined as the following infinite 



Mu,p) =Mu) = -i E (-l)V'""^'^'/^e(2-+i)i" = 2 X](-l)V'"^'^'^'sin(2n+ 1) 



sums i'2 : 












^2{u,p) = 


^2 




^3(m,P) = 


^3 




'diiu,p) = 


1^4 





u 



n=0 



(271+1)2/4 



cos(2n + l)u 



71= — OO 

oo 



71= — OO 
OO 



n=0 



1 + 2 Ep"'cos2nM 



(C.l) 



n=l 



COS 2nu. 



71=1 



By virtue of Jacobi's triple product identity |TP[, the theta functions admit a representation 
as infinite products 

OO 

{u) = 2 p^l^ sin M n (1 - V" cos 2u + p'") (1 - 

n=l 



i?2(m) = 2pi/^ cosu n (1 + V" COS 2u + - 

n=l 

oo 

|?3(m) = n(l + V""'cOS2li+p^"-2)(l-/n) (C.2) 

n=l 

oo 

n=l 

Another function which proves to be useful is 

oo oo 

n=— oo n=l 

Prom its definition it follows immediately that 

E{x,p) — E{px~^,p) — —xE{x~^,p) 

E{xp'',p) = (-a;)->-"("-')/'£;(x,p) 
E(p,p') = Q{p) (C.4) 

where the function Q is defined as 

oo 

Q{p)=Ui^-Pn- (C.5) 

n=l 

The ■j?-functions can be expressed in terms of the ^'-function as 

i)s{u,p)=E{-pc^'^,p^) (C.6) 

Likewise, the conjugate modulus transformation of the theta functions, which relates 
functions of nome p — exp(— e) to those of nome 

p' = e-'^'/^ (C.7) 
can be written as 

^2{u,p) = (^] e-" /^£;(p'e-2™/e,p'2) 

/^\ 1/2 

Mn,p) = (^j e-"'/^£;(-p'e-'""/Sp'') (C.8) 
= (^^y^'e-(«-/2)Ve^(_e-2W^^p'2)_ 

We note that these relations follow directly from the Poisson summation formula. 
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